APPLIED  PHYSICS  LABORATORY 

M>2!\GeDrR?a  A.e'Vi':  •  S".-'-r  .Sp-'fc  ■.Va,.<j')'! 

•  .  T  »•«'*.  w 


TRANSLATIONS 


CLB-3  T-649  2  December  1971 


OPTIMAL  PARAMETERS  OF  A  SATELLITE 
GRAVITATIONAL  STABILIZATION  SYSTEM  WITH  GYRO-DAMPING 


00 

'£> 

N 


V.  A.  Sarychev  and  K.  V.  Lukanin 


translated  by  L.  Hcltschlag  from 


Ordena  Lenina  Institut  Prikladnoi  Matematiki  Akademl i  Nauk  SSSR 
[the  Order-of-Lenin  Applied  Mathematics  Institute  of  the  Acad.  Sci.  USSR] 

Preprint  No.  46,  Moscow,  1971 


SUMMARY 


This  paper  investigates  the  transient  attitude  performance* 
of  a  gravity  gradient  satellite.  Damping  is  provided  by  a  pair  of  »wo- 
degree-of-freedom  gyroscopes.  First  the  dynamical  equations  of  motion 
are  presented  and  used  to  obtain  the  characteristic  equation.  The 
parameters  which  optimize  the  spacecraft's  transient  attitude  perfor¬ 
mance  are  determined  from  the  solution  of  the  characteristic  equation. 
The  results  are  then  specialized  to  the  case  of  an  axisymmetric  satel¬ 
lite  and  to  the  case  of  an  external  disturbance  due  to  orbital  eccen¬ 
tricity. 
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Introduction 

The  dynamics  of  a  satellite  gravity-gradient  stabilization  system  with 
gyro-damping  are  examined  in  this  paper.  The  gyro-damper  consists  of  a  pair  of 
two-degree -of- freedom  gyroscopes;  the  axes  of  rotation  of  these  gyroscopes  in  the 
satellite  equilibrium  position  are  oriented  symmetrically  relative  to  the  normal  tc 
the  orbital  plane.  The  natural  satellite  oscillations  bring  about  precession  of 
the  gyroscope  rotors  coupled  to  the  damping  device,  leading  to  dissipation  of  the 
energy  of  the  system. 

^  Translated  from  Ordena  Lenina  Inst.  Prik.  Mat.  Akad.  Nauk  SSSR  FThe  Order-of- 
Lenin  Appl.  Math.  Inst.,  Acad.  Gci.  USSR],  Preprint  No.  46,  Moscow  (1971). 
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References  [1-8]  have  been  devoted  to  the  study  of  such  systems.  Insofar 
as  the  range  of  problems  examined,  the  present  paper  is  closest  to  [1],  Thus,  with 
respect  to  determination  of  the  maximum  speed  of  response,  the  results  of  the  present 
paper  are  in  a  particular  case  a  refinement  of  the  results  of  [1]. 

Greatest  attention  is  paid  in  this  paper  to  determination  of  the  para¬ 
meters  with  which  maximum  system  response  time  is  achieved.  The  cases  of  an  axi- 
symmetric  satellite  and  of  a  satellite  with  arbitrary  distribution  of  mass  are  investi¬ 
gated.  In  the  study  the  values  of  the  optimal  parameters  were  found  by  the  random 
search  method  with  subsequent  refinement,  based  on  a  hypothesis  regarding  the  nature 
of  the  location  of  the  roots  of  the  characteristic  system  equation  at  the  point  of 
extremum.  In  the  neighborhood  of  the  optimal  parameters,  level  lines  of  the  degree 
of  system  stability  are  constructed  as  a  function  of  two  parameters,  the  remaining 
parameters  fixed.  The  dependence  cf  the  degree  of  system  stability  and  of  the  optimal 
parameters  on  the  amplitude  of  the  satellite  eccentricity  oscillations  is  constructed 
for  an  axisymmetric  satellite. 

1.  Equation  of  Motion 

The  configuration  of  a  satellite  gravitational  stabilization  system  with 
gyro-damping  is  described  in  detail  in  [1-3].  The  satellite  body  contains  two 
two-degree-of- freedom  gyroscopes,  which  are  used  to  provide  for  damping  of  the  natural 
oscillations  of  the  system  and  to  obtain  additional  restoring  moments. 

l^er  us  define  the  right  rectangular  coordinate  systems  necessary  for  the 
solution  of  the  problem. 

OXYZ  is  the  orbital  coordinate  system.  The  OZ-axis  is  directed  along  the 
radius-vector  connecting  the  centers  of  mass  of  the  earth  and  sate lllto -gyroscope 
system.  The  OX-*xis  coincide#  with  the  transversal  to  the  plane  of  orbit,  the  0Y- 
axis  with  the  normal. 


Oxyz  i#  a  coordinate  system  linked  to  the  satellite  body.  The  Ox-,  Oy-, 
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Oz -axes  are  the  principal  axes  of  inertia  of  the  satellite-gyroscope  system.  0  is 
the  center  of  mass  of  the  system. 

Qi  xyz  and  (fexyz  are  coordinate  systems  whose  axes  are  parallel  to  the  axes 
Oxyz;  <\  and  (fe  are  the  centers  of  mass  of  the  gyroscopes. 

The  axes  of  the  gyroscope  gimbal  suspensions,  a'b'  and  a"b"  {Fig.  1),  are 
parallel  to  each  other,  are  located  in  the  planes  xtfez  and  x(fez,  reap.,  and  form  the 
angle  Xq  with  the  axes  (fez  and  (fez.  The  axes  of  rotation  of  the  gyroscopes  In  the 
equilibrium  position  in  circular  orbit  form  the  angle  6Q  with  the  axes  (fey  and  (fey. 

The  location  of  the  satellite-linked  coordinate  system  relative  to  the 
orbital  coordinate  system  is  defined  by  means  of  the  angles  a  (pitch),  3  (yaw)  and 
Y  (roll)  (Fig.  2).  The  elements  of  the  direction  cosine  table 
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are  expressed  in  terns  of  these  angles  in  the  following  manner; 


am  -  cosu  -ffsji . 

Cfa  -  Sin  Si  f  -  Cl  :-  1  Sin/.'  easy , 

G,t  *  iu ■  t  c. ■:>  r »  CCStf  Sinji  Sr  n  jf , 

Qu  -  iVsn  CCSg  , 

Gu~  Cl  i  jSi  •>  i  »  Sin  jtSinfl  cq  >  J  , 

<7„ *  -v  Si  ;*JSinfi  $i.i  r  . 

(1) 

In  [2]  the  equations  of  motion  of  the  satellite-gyroscope  system  were 
derived  in  an  exact  formulation  without  resorting  to  the  precession  theory.  Assuming 
that  the  characteristic  rate  of  gyroscope  rotation  is  high  and,  consequently,  that 
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iaiifirtrtirViiii  gfiifiii<Mg~f£ft*t8Yffi' 


*od  the  dimensionless  time  T  *  oi0t.  Her 
«0  *  /fl/P3 »  where  £  is  the  grevitational 
constant  and  p  is  the  orbit  semi-major 
axis. 


Fig.  2.  Relationship  between  the 
trihedron  Oxy*  and  the 
orbital  coordinate  system. 


Then  the  equations  of  motion  of  a  satellite  with  gyroscopes  has  the 


form  [ 2] 


9Ap  =  H(pa 8,  *fi„ A J  >/„)]- 

“  9  /  $  *  ♦  H(  fa  *pn  j-  3(U  eeos  vf(l-  Qc )  an  an , 

<j  =  4  *fi'L  4 )  * pfect+WAn  * jL  )]  - 

-  ?fiAp*H(/in*/nj+  i(1*ecosvf(6A-Qe)an  a„  , 

6ei  =  H(/iu  S, */n <$z)*<i[6tp* H(p,i * pn )j  - 

-  pfe  *  H(pn  *  pn )]  *  3( h  ecos  vf(i-  Ba  )  a},  ait  , 

K  &  *-a •tS,*H$inS0  -  H(ppn  *qpu  *  ?jin  )  , 

Kt  Sz  =-Ai  Sz-HsinSe  -  H(ppn+<lfln  *  %fin  )  . 


Here  the  dot  denotes  differentiation  with  respect  to  dimensionless  time  t. 


fin  •  cos  X,  ccs(S,  *  6, ), 


fin‘cosX,cc$(-S\.  St). 


fin  --  -  cos  4  sinfS,  *  fi,)t  fil}  -  -  COS  Ac  sin(-  S,  *  St } 

I 

'l:  ~  S~n($t*  A,  I <  fill*  St/fl-  Or*  d.i  , 


fltl*ce$(S0'&). 

Pti r  ~sin  A,  cos  (6,  *  6,). 
fin  -  Sin  At  SinfS,  ■  c.)' 


fin 1  cos(-  <£  *6») . 

fin  •- sin  A,  cos  (-&»$*). 

fin  -  Sin  A,  sin  (-  A,  *  o,), 


(3) 


p,  q,  r  are  the  projections  of  the  absolute  angular  velocity  of  the  satellite  onto 
the  axes  of  the  coordinate  system  Oxyz,  and  65  are  the  angles  cf  deviation  cf 
the  gyroscope  rotational  axes  relative  to  their  equilibrium  positions  i-»  a  circular 
orbit,  v  is  the  true  anomaly  and  e  is  the  orbital  eccentricity. 

The  relationship  between  the  projections  of  the  absolute  angular  velocity 
of  the  satellite  (p,q,r)  and  the  angles  a,  g,  Y  is  given  by  the  formulas  [2] 


MalffiflMiifflfrT-tTrigft  in 
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Pm  *atl*  j*  Ca„  . 

9  «■  dan  *  fi  sir.#  *  valt , 
7  *  c(otK  -  ji  cos/f  *  vati  . 


where 


Vi(i*  ecos  \t  f  , 


(5) 


In  tht  circular* orbit  case  the  system  of  Eqs.(2j-(4)  has  the  solution 
1  =  0  *  Y  *  6i  *  62  "  p  *  -  *  0,  q  *  1,  which  corresponds  to  a  satellite  which  has 
no.  mo*-  ’d  with  respect  to  the  orbital  coordinate  system.  The  system  equations  of 
»:>ot<on  i  •  a  Ircular  orbit,  linearised  in  the  neighborhood  of  this  equilibrium  posi¬ 
tion,  ;i'  cribed  in  tht  following  manner,  after  introducing  the  notation  * 

H  sin  Oq,  Hg  *  H  cos  6g,  K  *  +  H2 : 

o(*3:Gt-et)o(  -  H (#,-&)*  0. 

2H,6t  ♦  *,(6,-61)  *■  M(d,-St)  -  0 ,  (6) 

Qck  *■[£ Ht -,«?«*  0C  -  t)]j  *[2Ht  * 

♦  ♦  S.  ' r.r.L*  H.(6,*6}  )  ca  A 0  -  0  . 

+[2Ht<*(/-9')}f- 

~Hl(6,*6t)co&Ae  +  Ht(6,*6t)sbi  A„  =  0. 

'•2HtfisLnAt  +  2Xxfic*sA,*  co$Xe  - 

*2Ht  A,  *■  kj<5,‘8t)+  v(S,  +  ot)  ■*  0.  (7) 

2.  Method  of  Determining  the  Optimal  Parameters 

me  of  the  criteria  for  the  quality  of  stabilization  systems  is  the  time 
required  for  transient  processes  to  take  place  in  them.  J.u  the  present  paper  it 
was  assumed  Chet  the  transient  response  time  is  defined  by  the  variable  |,  which  is 
the  distance  between  khe  rightmost  root  of  the  characteristic  equation  of  the  system 


6* 


(6)-(7)  and  the  imaginary  exit.  The  characteristic  equation  of  system  (6)-(7)  is 


written  in  the  form  of  the  product  of  the  third-  end  fifth-degree  equations 


■a,X.(a±-a,*£  ).  0. 


Sin  A,  COS* X„ 

a'--er  +  ~  • 


.  '*(i-BtbfWz-(9»+Qt •  it 

eA  o. 


„  . i(8* * et-ihZHz •  (4- Sic )sinA,  +  (t-20i)cos*L 

eA0t 


a  -  Sin* A* *  4(f-%  jcn*A, 


The  variable  »  depends  on  the  seven  dimensionless  parameters  9A,  -c,  H, 

,  Kq  ,  5q»  X.q.  The  problem  consists  in  determining  the  values  of  these  parameters 


at  which  Z  reaches  maximum  value  ®  .  These  parameter  values  will  be  said  to  be  optimal. 


the  present  paper  we  are  studying  net  only  a  satellite  with  arbitrary  moments  of 


inertia,  but  also  the  case  of  an  ay.isyametric  satellite  (a  *  1). 


The  investigation  was  carried  nut  in  two  steps.  The  first  step  consisted 


in  determining  *  and  the  corresponding  optimal  parameters  by  the  random  search 


method.  In  the  second  step,  the  parameter?  we  found  were  determined  more  exactly 


on  the  basis  of  a  hypothesis  as  to  the  form  of  the  roots  of  the  characteristic  equa< 


tion  at  the  point  of  extremum. 


The  essence  of  the  random  search  method  is  as  follows.  We  begin  with  a 


point  Pq  in  the  parameter  apace  of  r.he  system.  Frost  this  point  we  take  the  tandow 
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*JL*- i I < ^ 1W^ J PJ  l j . '.U'  -...'.I.- - •  — rr— — 7 - — - — .  ~ 

step  Apj  .  If  §(pQ  +  Apj  )  >  '(pg),  then  the  point  pQ  +  Apj.  is  taken  as  the  initial 
point  for  the  next  random  step.  If  ?(pg  +  Api  )  £  J(pg),  then  we  take  the  random 
step  Ap2  from  the  point  pg.  The  point  pg  was  considered  to  be  the  extremal  point 
if  the  variable  §  did  not  increase  in  ra  random  steps  from  this  point. 

The  value  of  5  was  calculated  with  the  prescribed  accuracy  by  the  method 
of  successive  approximation.  The  zero  approximation,  the  interval  containing  the 
rightmost  root  of  Eqs.  (8)-(9),  was  determined  by  the  formula  proposed  in  [l].  The 
subsequent  refinements  were  carried  out  by  the  Routh  scheme.  Routh's  scheme  is 
constructed  for  equations  obtained  from  (8)-(9)  by  right-displacement  of  all  of 
their  roots  by  a  value  equal  to  the  distance  between  the  centerpoint  of  the  interval 
containing  %  and  the  imaginary  axis.  A  detailed  description  of  the  process  of  cal¬ 
culating  "  is  given  in  [ l]  . 

In  the  process  of  determining  by  the  random  search  method,  it  was 
discovered  that  if  the  values  of  "m  obtained  at  the  end  of  optimization  from  various 
initial  points  differ  by  several  percent,  the  scatter  over  the  parameters  themselves 
curns  out  to  be  considerably  greater. 

The  difficulties  involved  in  precise  determination  of  the  optimal  para¬ 
meters  are  due  primarily  to  the  great  size  of  the  parameter  space  and  to  the  high 
degree  of  the  characteristic  equation.  Moreover,  as  the  investigation  showed,  s 
not  insignificant  part  is  played  by  the  nature  of  the  extremum  itself,  near  which 
the  derivatives  of  *  with  respect  to  several  directions  in  parameter  space  differ 
substantially  in  value.  This  is  explained  by  the  fact  that  the  real  parts  of  the 
roots  of  the  characteristic  equation  become  equal  near  the  extremum;  as  a  conse¬ 
quence,  the  function  Z  terns  out  to  he  non* analytic  at  the  extremal  point. 

3 .  Optimal  Parameters  of  a  Satellite  with  Arbitrary  Moments  of  Inertia 

The  satellite  parameters  were  optimized  by  the  random  search  method  for 
fixed  values  of  the  angle  \g.  As  mentioned  above,  the  values  of  obtained  at 
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the  end  of  optimization  from  various  initial  points  differ  by  several  percent,  while 
the  scatter  over  the  parrmete ra  turns  out  to  be  much  greater.  Given  in  the  table 
are  the  maximum  values  obtained  for  for  each  of  the  fixed  values  of  Xq,  as  well 
as  the  corresponding  roots  of  characteristic  equations  (8) -(9). 
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In  this  table  yx  ,  ys,  y3  are  the  roots  of  Eq.  (8),  xl ,  Xa,  jfe  ,  x* ,  Xs  are  the  roots 
of  Eq.  (9). 

from  the  table  it  can  be  seen  that  the  real  parts  of  the  roots  y2 ,  ya  , 

X2 ,  Jfc  ,  x*  ,  x*  become  approximately  equal  in  the  neighborhood  of  the  extremum  *, 
and  with  \q  *  O'5,  30° ,  and  6<f'  the  root  also  becomes  equal  to  the  real  parts  of 
these  roots.  The  investigation  was  begun  with  the  most  complex  case  X.-,  *  90°. 

J 

Let  us  advance  the  following  hypothesis  as  to  the  nature  of  the  roots-  at 
the  point  of  the  real  extremum  -  with  \q  *  90°  the  roots  of  the  characteristic 
equations  (8)- (9)  have  the  form 

&  x'0'' »  H*. i <«» ,  JTi *  - <**  -li.,  «=- eft  (ut 

iu)j  ,  q (,*</,  oi  *  </. 

(ID 

The  variables  d,  ,  da,  ju^  ,  ufc,  are  subject  to  determination.  Naturally,  of 
greatest  Interest  is  the  parameter  d,  since,  by  construction,  d  ■ 

The  characteristic  equations  with  roots  (11)  have  the  form 


rAftiuwtrtjiittWri-i 


I 


(12) 


AJ*  (Zd*d,)/*[(di  <4  )  +2dd,]k  *  df(^  d?)  -  0  . 

A  *  (<id*  dt  )A  ** [(ul*  co?J  *■  4ddt  *  6d*JAs* 

*[(2d*dtXol<rut)  *  6d1cti+<rct]A*+  [tf(c(*2dtX<4+<>>!)+ 

■no*,u£~4afa'1  *d*]A *  dt [d*(a£*uf ) *■  c£aif*d*] *  0.  (1j) 

Let  us  impose  constraints  on  the  parameters  of  Eqs.  (8) -(9)  so  that  the  roots  of 
these  equations  will  always  have  the  form  of  (11).  ?o  do  sot  we  equate  the  coeffi¬ 
cients  of  Eqs.  (8)-(9)  to  the  corresponding  coefficients  of  Eqs.  (12)-(13),  and  we 
also  rewrite  the  last  two  inequalities  of  (11).  Then, 
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a.  4 

K,  **. 

3(Ga-Q')  =  c^c,'%2g3-.  . 

-J  3(9^6.) *ct.f4*** 

*  4d*a- 

K  if,  v  •  ■ 

a.  -  fa:** cf. } -  'fci.\  -  {-  a* . 
a.~*at ±2*  -- (2d*t f. jfcjf *a)J)*Gdif, * 4# 3 , 
as:'d(d-2 c\ Xu, -of,)*  cjjol* 4c'cft  *  d*, 
djd%c%<4 ) > c4<£* d*J . 
d,*d  .  dx  *  d . 


From  (14)  we  get  the  following  expressions: 


<S,*3(GA- %)-</*•  2dd„ 

*  d,(d**  uf)  , 

3iGa-8.) 


~2d*d,-£ 


u\*a)j*err- 4dda  ~  6d* . 

•  ard(d*  2dtX*{*c£)-«dadj  "  d ' 

o5*c(, 


'  .?<*  * 
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(14) 


(15) 

(16) 


(a'-!kh‘  <”> 

Here, 

2t*(4d*dt  )%r  +dt[dt(<4*i-il )*<£<£  *  d*J, 

V -(«d*dt)%t*(2d*dk)(ut*vl)+  Sdcft  +*fd\ 

T\  9  /tl  C  9(dk~  Sc)z< 

-  (U) 

If  relations  (15}-(18)  hold,  the  roots  of  the  equation?  being  investigated,  (8)-(9), 
can  be  represented  in  the  form  of  (11). 

By  imposing  constraints  (15)-(18),  a  set  containing,  by  hypothesis,  the 
optimal  parameters  has  been  singled  out  in  the  space  of  the  parameters  9A»  Gc»  H, 

Kj ,  Kq ,  6g.  These  optimal  parameters  correspond  to  the  maximum  value  of  d  at  which 
all  the  variables  in  (15) — ( 18)  satisfy  the  physical  meaning.  Therefore,  in  addition 
to  the  purely  formal  constraints  specified  by  (15) - (18) »  it  is  also  necessary  to 
require  fulfillment  or  the  conditions  uf  i  0,  u§  £  0,  u£  £  0,  D  £  0,  £  0,  1%  £  0. 

It  is  also  necessary  to  remember  that  the  sum  A  two  moments  of  inertia 
is  always  greater  than  the  third,  and  therefore  9A  +  GC  *  ®*nce  Ki  >  0  in  the 

case  of  asymptotic  stability,  the  conditions  £  0,  H|  £  0  reduce  to  the  requirement 

2H?  2H|  _  , 

that  the  quantities  -rr-  be  non-negative.  Moreover,  £  0,  ug  £  0  and  are 

real.  Consequently, 

6^  +  6^  *0,  *  (?,  (<*>•*  (ty)  -  9cj*u*  *■  O. 


Thus,  in  order  for  the  solution  of  the  characteristic  equations  to  exist  in  form  (11), 
it  is  necessary  to  require  fulfillment  not  only  of  relations  (15)-(18),  but  also 
of  the  inequalities 

mI<*\  *  0,  (a}\ *u[f-  a  O  ,  Bk+Qc  +  »,  (19) 


At  a  fixed  value  of  the  angle  \q,  expressions  (10)  for  a^,  a1  ,  ^ ,  as,  aA 

depend  on  the  variables  0A,  9C,  Hg ;  rx  and  r2  in  (18)  aepsnd  on  the  variables  9A, 

°q>  Ha,  d2 ,  d.  Therefore,  relation  v17)  can  be  used  to  define  da  m  termo  of  9A, 

^c>  Ife ,  d.  Considering,  moreover,  that  dx  is  calculated  by  formula  (16),  we  get  that 

it  is  possible  to  choose  four  parameters  arbitrarily  in  relations  (15) - (18) :  aA, 

Ife ,  d.  These  parameters  must  be  chosen  such  that,  by  satisfying  inequalities  (19), 

we  get  the  maximum  value  of  d,  equal  in  this  case  to  ?  . 

m 

The  domains  in  which  inequalities  (19)  and  the  inequalities  dx  s  d,  da  >  d 
are  fulfilled  were  constructed  in  the  plane  of  the  parameters  9A,  9  (Fig.  3)  for 
fixed  values  of  Hq .  The  variable  d  increased  steadily  until  the  corresponding  do¬ 
main  contracted  to  a  point.  Plotted  in  Fig.  3  are  the  boundaries  of  the  domains 
for  d  *  0.5,  0.52,  0.525.  As  can  Le  seen  from  the  figure,  the  domains  of  existence 
of  a  solution  in  form  (11)  for  any  H2  are  determined  by  the  curves(ii|  4-  rji£);'  - 
4  *  0  (curve  I),  D  >  0  (curve  II),  dg  -  U  (curve  III)  and  the  straight  line 


Fig.  3.  Domains  in  which  the  inequalities  ar*  fulfilled  for  fixed  values  of  Ife  . 


JA  +  9 £  ■  1.  The  sign  in  front  of  the  root  in  formula  (16)  for  d,  has  no  effect  on 
the  form  of  the  domains;  it  merely  brings  about  a  double  value  in  the  definition  of 
9^  .  The  maximum  value  d  as  0.525  is  reached  when  Hg  «  0.19. 

A  more  important  result,  however,  is  the  fact  that  for  any  1^  there  is  a 
curve  dg  “  d  (curve  III)  among  the  curves  bounding  the  domains  of  existence  ot  roots 
in  form  (11).  Consequently,  the  equality  dg  *  d  will  hold  at  the  real  extremal  point 
For  values  of  >.q  «  0°,  30°,  60°  this  was  evident  from  the  results  obtained  by  the 
random  search  method.  Thu3,  for  any  \q,  the  inequality  d2  &  d  in  (11)  can  be  re¬ 
placed  by  the  equality  d2  *  d.  This  makes  it  possible  to  choose  a  total  of  three 
parameters  arbitrarily  in  formulas  (15)-(18):  9.,  9  ,  d. 

A  t> 

In  (15)-(18)  we  set  da  =  d.  Then  relation  (17) t  which  contains  9A,  9C» 
d,  Hg,  can  be  used  to  determine  Hg  for  given  values  of  0^,,  d.  From  (15)  we 
find  in  addition  that 

,  / 

A,  _ 

*» ♦  T/j  j  6t*azccc:(yj  (20) 

Formulas  (20)  associate  definite  values  of  the  parameters  H,  6q,  1^  ,  Kg  with  any 
point  in  the  domain  of  existence  of  the  solution  in  the  plane  of  the  independent 
parameters  9^,  9^,;  at  the  same  time,  because  of  the  two  values  in  formula  (16)  the 
variables  H  and  trj  will  also  take  on  two  values.  This  double-value  situation  disap¬ 
pears  on  the  curve  D  *  Q.  Formulas  (20)  make  it  possible  to  map  the  boundaries  of 
the  domains  of  existence  of  the  solution  from  the  plane  of  the  independent  para¬ 
meters  9a0£  onto  the  plane  of  one  dependent  and  one  independent  parameter  ,  e.g.s 
9.-  Such  a  mapping  ia  convenient,  in  that  it  gives  us  some  idea  of  the  nature  of 
tiie  function  *  in  the  region  of  its  extreraut*. 

The  results  obtained  for  values  of  XQ  *  OP,  30°,  60°,  90°  are  presented 
in  Figs.  4-13.  As  can  be  seen  from  these  figures,  the  domains  of  existence  of 
solutions  of  form  (11)  with  d~  »  d  are  determined  by  two  curves: 


u‘\  ./ 

\  *• ' 
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(u£  +  U§)3  -  4  u|(jJ  «  0  (curve  1),  0  *  0  (II)  and,  beginning  with  a 
certain  value  of  \q,  the  line  0^  +  •  1  (III).  Because  dx  is  two-valued  in  its 

definition,  curves  I  and  III  of  the  plane  of  9^,0^  are  mapped  into  the  curves  l+, 
l",  III+,  III~  on  the  planes  of  6>f  H  and  9^,  5^.  D  *  0  on  curve  II,  and  therefore 
the  definition  of  d-  is  not  two-vai»ed.  The  cross  in  all  the  figures  denotes  points 
at  which  the  maximum  value  d  ■  is  reached.  At  the  optimal  point,  D  *  0,  and 
therefore  H  and  6q  are  uniquely  determined.  From  Che  nature  of  the  curves  d  •  const 
it  can  be  seen  that  in  the  region  of  the  extremum  the  function  ?  is  less  sloping 
when  \q  «  0.  With  increasing  \q  the  distance  between  the  curves  d  ■  const  decreases. 
The  considerable  elongation  of  the  figures  described  by  the  curves  d  *  const  con¬ 
firms  the  fact  found  by  the  random  search  that  the  derivatives  of  £  with  respect  to 
various  directions  in  parameter  space  differ  substantially. 


Fig.  14.  The  variable  and  the  shape  of  a  satellite  in  the 
c  hiatal  configuration  as  a  function  of  Xq. 
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Fig.  15.  Dependence  of  ~A>  «c,  H,  Kj ,  ,  5q  on  XQ. 

Shown  in  Fig.  14-15  is  the  dependence  of  and  of  the  optimal  parameters 
§A»  9r»  H,  Ki  ,  !fe,  6q  on  Xq  >  1  when  Xq  <  27°,  and  the  optimal  satellite  con¬ 

figuration  is  gravitationally  unstable,  since  the  condition  1  >  9.  >  dr  should  be 
fulfilled  for  stability.  Stability  is  achieved  by  the  kinetic  moments  of  the  gyro¬ 
scopes.  The  corner  points  on  the  graphs  for  \n  m  35°  are  explained  by  the  fact  that 


1.  ]‘.n 


when  Xq  a  35°  the  optimal  parameters  are  found  on  the  boundary  9^  +  9C  * 

this  case  the  satellite  is  a  gravitationally  stable  plate.  We  note  that  when  XQ  > 

31°,  the  optimal  value  of  the  elasticity  coefficient  Kg  is  less  than  zero.  The 

maximum  value  «  0.545  is  reached  when  \Q  -  90°,  $A  -  0.857,  a  *  0.143,  iq  * 

0.238,  Kg"  -0.078,  H  *  0.476,  5q  *  65.78°.  A  variant  of  the  configuration  for  * 

9CP  was  also  considered  in  r 1] ,  where  ?  *  0.48  was  obtained. 

■ra 

The  optimal  parameters  9^,  9C»H,  iq  ,  Kg,  5q  were  sought  on  a  set  yielding 
roots  of  form  (11)  for  da  «  d.  in  order  to  prove  that  the  parameter  values  we 
found  do  indeed  yield  a  local  maximum  of  it  is  necessary  to  study  the  neighbor¬ 
hood  of  the  point  of  extremum  for  arbitrary  parameter  values.  The  following  con¬ 
sideration  turns  out  to  be  convenient. 

We  fix  four  of  the  six  optimal  parameters  and  in  the  plane  of  the  remaining 
two  we  construct  a  region  within  which  the  degree  of  stability  *  is  greater  than  a 

given  value  of  *Q.  To  do  so,  we  make  the  exchange  \  ■  »Q  +  o  in  Eqs.  (8)-(9).  The 

region  to  be  sought  is  one  in  which  all  the  Hurwttz  determinants  of  the  character¬ 
istic  equations  obtained  with  respect  to  o  are  positive.  By  assigning  several 
values  to  we  get  a  number  of  regions  imbedded  in  each  other  with  boundaries 
yielding  level  lines  of  “  as  functions  of  the  two  chosen  parameters  for  optimal 
values  of  the  remaining  parameters.  The  closure  of  the  level  lines  of  -  will  be  a 
proof  that  the  parameter  values  we  found  do  indeed  yield  a  local  maximum  for  f. 

Constructed  in  Figs.  16-19  are  level  lines  corresponding  fc  values  of 
%  *  0.2,  0.3,  0.4  for  all  parameter  combinations  with  *  90°.  The  cross  denotes 
the  point  of  maximum.  As  can  be  seen  from  the  figures,  all  the  level  lines  of  * 
are  closed  and  encompass  the  point  of  maximum. 

Figs.  16-19  give  an  idea  of  the  behavior  of  the  function  *  in  the  region 
of  the  maximum.  They  may  turn  out  to  be  useful  in  a  structural  realization  of  the 
stabilization  system,  when  it  is  difficult  to  get  precise  optimal  parameter  values. 
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4.  Optical  Parameters  of  an  Axisymmetric  Satellite 
An  axisymmetric  satellite  is  defined  by  the  condition  »  1.  A  study 
made  by  the  random  search  method  has  shown  that  in  this  case  the  equality  d2  *  d 
does  not  hold  for  all  values  of  Xq.  In  particular,  for  X  *  90°,  this  is  clearly 
evidenC  in  Fig.  3.  In  fact,  if  for  any  Hg  we  decrease  d,  then  for  a  certain  valte 
of  d  the  straight  line  3^  *  1  intersects  curves  I  and  II,  whereas  curve  III  (da  *  d) 
will  remain  on  the  left.  This  reasoning  forces  us  to  seek  the  optiasl  parameters 
on  a  set  specified  by  conditions  (11).  In  (15)-(18)  we  set  9  »  I.  Relation  (17), 

A 

which  contains  9g»  Ha,  d,  da,  can  be  used  once  again  to  determine  da.  Then,  by 
formulas  (20),  w*  find  all  the  remaining  parameters.  Although  in  this  case  the 
independent  parameters  are  3g  sod  Ha ,  it  is  only  natural  not  to  construct  regions 
defined  in  the  plane  of  these  parameters  by  the  Inequalities  (19)  and  dj^  2;  d, 
da  a  d,  but  to  proceed  immediately  to  construct,  by  formulas  (20),  a  mapping  of  the 
boundaries  of  these  regions  onto  the  planes  of  one  of  the  dependent  parameters  H, 

Ki ,  Kg ,  6q  and  of  the  independent  parameter  @c- 

The  results  obtained  for  values  of  Xq  “  0°,  30°,  60°  and  90°  are  shown 
in  Figs.  20-24.  When  9^  «  1,  the  domains  of  existence  of  the  solution  in  form  (11) 
are  defined  by  three  curves:  (u|  +  )a  -4  u|uf  *  0  (curve  I),  D  ■  0  (curve  II), 
and  d2  *  d  (curve  III).  The  nature  of  the  domains  changes  substantially  for  various 
Xq.  If,  however,  for  Xq  «  0°,  30°  it  still  turne  out  that  da  *  d  ^  at  the  point 
of  optimum,  then  for  Xq  *  60°  the  curve  da  ■  d  occurs  only  when  d  *  0.42,  while 
the  point  of  optimum  is  defined  by  curves  I  and  II,  and  d2  >  d  *  at  this  point. 

As  in  the  case  of  arbitrary  values  of  9^,  when  9^  “  1  the  figures  described  by  the 
curves  d  •  const  are  seen  to  be  quite  elongated. 

The  dependence  of  and  of  the  optimal  parameters  9C,  H,  1^,  Kg,  on 

Xq  is  plotted  in  Fig.  25.  The  maximum  value  5^  *  0.483  is  reached  when  Xq  *  9G° , 

6C  *  0.172,  ^  «  0.501,  ife  *  0.018,  H  -  0.703,  5Q  -  64.6°.  The  level  lines  of  - 
for  this  case  are  shown  in  Figs.  26-29. 
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5.  Eccentricity  Oscillations 

In  an  elliptical  orbit,  due  to  the  non-uniformity  of  rotation  of  the 
radius-vector  connecting  the  centers  of  mass  of  the  earth  and  the  satellite,  the 
stable  equilibrium  position  being  studied  goes  over  into  periodic  oscillations,  called 
eccentricity  oscillations.  The  period  of  the  eccentricity  oscillations  is  equal  to 
the  p< riod  of  revolution  of  the  satellite  center  of  mass  along  the  elliptical  orbit; 
the  oscillation  amplitude  is  proportional  to  the  magnitude  of  the  orbital  eccen¬ 
tricity. 

The  eccentricity  oscillations  are  a  particular  solution  £5  *  y  *  0,  i- 
^2  s  0,  a  ■  a(t),  5i  -  bz  *  f(t)  of  system  (2).  This  particular  solution  must  be 
found  in  the  form  of  a  series  in  powers  of  the  ec;  entricity  e.  Limiting  ourselves 
to  first  terms  in  the  expansion  in  e,  we  find  that  the  amplitude  Rj  of  the  satellite 
eccentricity  oscillations  is  determined  by  the  formula  [1,8] 
yp*  /  *  *sin*6,]Z 

'  *  e  «?[$(£* * Hcos ScfifcO'i-  *  (21) 

The  transient  response  of  the  satellite-gyroscope  system  in  the  case  of 
a  circular  orbit  is  constructed  in  Fig.  30.  The  satellite  is  axisymmetric  (9^  *  1) 
with  parameters  optimal  when  \q  ■  90°.  The  initial  values  of  the  angles  and  velo- 
cities  were  chosen  a  follows:  a  *  0  *  y  *  10°,  »  62  *  p  *  r  =  0,  q  *  1. 

From  the  figure  it  can  he  seen  that  in  two  satellite  orbital  revolutions  the  oscil¬ 
lation  amplitude  with  respect  to  ell  angles  becomes  less  than  one  degree. 

The  transient  response  of  the  satellite  under  the  same  initial  conditions 
In  the  case  of  an  elliptical  orbit  with  eccentricity  e  *  0.02  is  given  in  Fig.  31. 

The  initial  value  of  the  true  anomaly  v  was  chocen  tc  be  zero.  Steady^state  periodic 
motion  in  pitch  with  an  amplitude  of  3.6°  can  be  seen  in  the  figure.  This  value 
is  i.i  good  agreement  with  formula  (21),  which  yields  R  «  —  *  3.06,  from  which  we 
get  Rj  ■  R  e  ■  0.0612  rad  *  3,5°  for  e  »  0.02. 


. .  . . . . 


We  select  Che  parameters  of  Che  axisymnecric  satellite  being  considered 


such  that  maximum  E  is  obtained  for  a  fixed  value  of  R  *  ~  .  The  results  obtained 
5  e 

by  the  random  search  method  are  given  in  Figs.  32-33.  From  the  figures  it  can  be 

seen  tnat  the  value  of  R  can  be  reduced  from  three  units  to  two  without  decreasing 

the  degree  of  stability  appreciably.  When  R  *  1,  the  degree  of  stability  *  *  0; 

here  it  is  easy  to  show  that  this  will  remain  valid  for  a  satellite  with  any  desired 

value  of  \0  and  arbitrary  distribution  of  mass.  In  fact*  in  seeking  minimum  R  by 

formula  (21) ,  we  get  that  the  minimum  value  R  *  1  is  reached  when 

H' Sin' So  n 

X,  Z  X,  *'  * 

In  this  case,  on  the  basis  of  the  Routh-Hurwitz  criterion,  for  characteristic  equa¬ 
tion  (8)  we  find  that  ^  *  0.  The  point  R  ®  3.06  corresponds  to  the  optimal  satellite 
parameters  obtained  in  section  4.  The  break  in  the  graphs  for  R  *  3.06  is  explained 
by  the  differing  nature  of  the  roots  at  which  the  extremum  »  is  realized,  to  the 


right  and  left  of  this  point. 


Conclusion 


The  degree  of  system  stability  Z  is  not  an  analytic  function  at  the  point 
of  maximum.  This  is  due  to  the  fact  that  the  real  parts  of  the  roots  of  the  char¬ 
acteristic  equation  being  examined  become  equal  at  the  point  of  maximum  |.  Even 
though  it  has  negative  aspects,  this  fact  permits  us  to  make  the  optimal  parameter 
values  obtained  by  various  rough  numerical  methods  more  precise  by  introducing 
relations  between  the  parameters  obtainable  by  equating  the  real  parts  of  the  roots 
of  the  characteristic  equation  being  studied. 

As  applied  to  the  satellite  configuration  with  gyro-damping  being  considered 
here,  the  following  results  were  obtained  by  this  method.  The  maximum  value  of  the 
degree  of  system  stability  *  0.i>45  is  reached  when  \q  *  9v~ ,  which  corresponds 
to  ’©cation  of  the  gyroscooe  global  axes  along  the  Ox-axis;  in  the  optimal  configu¬ 
ration,  the  satellite  is  a  gravitationally  stable  plate.  When  \Q  <  Zf' ,  the  optimal 
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satellite  configuration  is  gravitationally  unstable,  and  system  stability  is  achieved 


by  using  the  kinetic  moments  of  the  gyroscopes.  When  Xq  >  31°,  the  elasticity  co¬ 


efficient  Ka  <  0.  The  maximum  degree  of  stability  of  an  axisymmetric  satellite  is 


also  reached  when  \n  ■  90°.  In  this  case  C  ■  0.483. 

u  wtn 


The  level  lines  of  5,  plotted  in  two  parameters  in  Figs.  16-19  £fu  an 


arbitrary  satellite  and  in  Figs.  26-29  for  an  axisymmetric  satellite,  allow  repre¬ 


sentation  of  the  function  §  in  the  region  of  the  maximum.  These  figures  may  prove 


to  be  useful  in  the  design  stage  of  the  system,  when  it  is  usually  difficult  to 


get  exact  parameter  values. 


The  study  we  made  of  the  dependence  of  ^  on  the  magnitude  of  the  amplitude 


of  the  eccentricity  oscillations  Rx  for  an  axisymmetric  satellite  yielded  the  follow¬ 


ing  results  for  \q  *  90°.  With  optimal  satellite  parameter  values, »  3.06e.  It 


is  possible  to  specify  parameters  for  which  Rx  *  2e,  and  the  degree  of  stability 


does  not  change  appreciably.  When  Rx  -  e  we  find  that  §m  “  and  it  -s  not  possible 


to  get  Rx  <  e.  This  is  also  trua  for  a  satellite  with  arbitrary  mass  distribution 


and  any 
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